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Abstract 

The critical behavior of two-dimensional stochastic lattice gas mod- 
els with C^y symmetry is analyzed. We study the cumulants of the 
order parameter for the three state (equilibrium) Potts model and for 
two irreversible models whose dynamic rules are invariant under the 
symmetry operations of the point group C3 V . By means of extensive 
numerical analysis of the phase transition we show that irreversibility 
does not affect the critical behavior of the systems. In particular we 
find that the Binder reduced fourth order cumulant takes a universal 
value U* which is the same for the three state Potts model and for 
the irreversible models. The same universal behavior is observed for 
the reduced third-order cumulant. 
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1 Introduction 



The critical behavior of nonequilibrium systems has been amply studied in 
the last years |], ^ ^ [|. These studies consider stochastic lattice models 
and probabilistic cellular automata that evolves in time according to an ir- 
reversible dynamics, that is, a dynamics that lacks detailed balance. An 
important aspect to be considered is the role of symmetry. Distinct systems 
with the same symmetry are expected to have similar critical behavior. The 
symmetry is to be found in the Hamiltonian for reversible systems and in the 
dynamics for the irreversible systems. Among the irreversible models, there 
are models that have a reversible counterpart with the same symmetries. In 
this context it has been established the following statement: models with 
up-down symmetry, similar to the Ising model, and defined on the same lat- 
tice, reversible or not, are in the same universality class 0. This has been 
verified numerically for a large number of models || 0, |], H [10], [Q[ We 
remark that the same universal behavior is also observed for dynamic phase 



transitions in Ising models in oscillating fields [13, 14, 15 



Recently |TB|, [T7|, [Tl|, we have argued that this statement can be extended 
to models with other symmetries. In fact, in these works, we considered prob- 
abilistic cellular automata with dynamics that possess C^ v symmetry and ver- 
ified that the values of the critical exponents, both static and dynamic, for 
irreversible systems are the same of the equilibrium three-state Potts. That 
is, irreversibility is irrelevant regarding the values of the critical exponents in 
systems with the symmetry of the Potts models. In this article we complete 
this analyzes by performing a systematic study of the cumulants at the crit- 
ical point. We consider stochastic lattice gas models with C% v symmetries: 
the three state Potts model ]nj and two irreversible models, and focus our 
attention in the determination of the third and fourth order cumulants of 
the order parameter. We present a Monte Carlo study of these properties for 
models defined in a regular square lattice. Our results show that the value of 
the cumulants, at the critical point, obtained for the present nonequilibrium 
models, and the ones associated to the two-dimensional three-state Potts 
(equilibrium) model are the same within numerical errors. 

The paper is organized as follows. In Section 2 we present the models 
to be studied. In Section 3, it is defined the quantities of interest in the 
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study of the phase transition and their scaling properties. The values of 
the cumulants in the limit of infinite temperature and zero temperature are 
discussed in Section 4. Sections 5 and 6 show numerical calculations and the 
conclusions remarks. 

2 Models 

Consider a regular lattice of iV sites in which each site can be in one of three 
states. At each site we attach a stochastic variable <7j that takes the values 1, 
2 and 3. The state of the system can be represented by a = (a±, Oi, • cr^). 
The time evolution equation for P(a,t), the probability of state a at time t, 
is given by the master equation 

j t P(a, t) = ^{W(a | a')P(a', t) - W(a' | a)P(a, t)} (1) 

where the sum is over the 3^ configurations of the system. W(a \ a 1 ) is the 
transition rate from a state a 1 to state a, given that at the previous time 
step the system was in state a'. We will consider dynamics where a and a' 
can differ only by one site which we call a one site dynamics. In this case we 
have 

Y f P ^ = E X>,(OiV Q , t) - Wi (a)P(a, t)} (2) 

where Wi(a) is the transition probability per site and 

a ia = ((7i,(72, ...,<Ji + a, a N ) (3) 

with a = 1, 2. The states are defined modulo 3. 

The transition probability Wia \ a') is invariant under certain symmetry 
operations, that is, W(Ra | Ra') — W(a \ a') where R is a symmetry 
operation. For the present models the symmetry operations R are those that 
act on all sites transforming each of them in the same manner. One of the 
symmetry operations is the rotation operation 1 — > 2, 2 — > 3, and 3 — > 1. 
Another is the operator 2^3 with state 1 fixed. If the three states are placed 
on the vertices of a equilateral triangle they correspond, respectively, to a 
rotation operation by 120 degrees and a specular operation. These symmetry 
operations define then the point group C 3v . 
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2.1 Equilibrium model and Metropolis prescription 

The Hamiltonian of the 3-state Potts model is 

# = £- J %i^-) (4) 

ft) 

where <7j = l,2,3,J>0is the interacting strength between the nearest 
neighbor sites i and j and 5 is the Kronecker delta. 

To simulate the model we associate to it a dynamics. We consider an one- 
site dynamics as described by the master equation (^). In the case the model 
to be analyzed is an equilibrium model it is necessary to use a transition 
probability Wi(o~) that satisfy detailed balance condition. That is, in the 
stationary state we must have 

w t (a ia )P(a* a ,t) = Wi (a)P(cT,t). (5) 

This dynamics can be defined by using the Metropolis prescription. We 
choose a site i and a state a and calculate 

Wi(a) = min{l, exp(-0AH)} (6) 

where AH = H(a) — H(a ia ) is the difference between the energy of the 
state a and the energy of the state a ta . The parameter (3 is associated 
to the inverse of the temperature T. Numerically we studied the critical 
point associated to the model by performing Monte Carlo simulations. This 
procedure is described in Section 6. 

2.2 Nonequilibrium models 

The nonequilibrium models are defined as follows. For the case of a square 
lattice we are denoting the transition probability w (a) by w(a \ai, <t 2 , 03, cr 4 ), 
where the sites 1, 2, 3, and 4 are the first neighbors of site 0. 

Symmetric stochastic lattice gas model 

(a) If in the neighborhood of a given site there is a majority of sites which 
are in one state then, independently of the state of the site, it changes to the 
state of the majority with probability p. It changes to one of the two other 
states with probability (1 —p)/2. 
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(b) If no state is in majority then the site assume either state with equal 
probability. 

According to the local rules of the model we have 

w(l | 1111) = w{l | 1112) = w(l | 1113) = w(l | 1123) = p 

and 

w(l | 1122) = w(l | 1133) = w(l | 2233) = 1/3 (7) 

The other rules are obtained by permutation of the neighboring sites and by 
cyclic permutation of the states. 

Majority stochastic lattice gas model 

The model consists of a stochastic lattice gas model where the site tran- 
sition probabilities follows similar rules of the majority vote model ||. The 
chosen site adopts the same value of the majority of the nearest neighbor 
sites with probability p. It adopts the state of the other states with proba- 
bility q/2 = (1 —p)/2. If there is an equal number of nearest neighbors sites 
in the same state then the chosen site adopts each state with probability p/2 
or it assumes the other state with probability q. That is: 

w(l | 1111) = w(l | 1112) = w(l | 1113) = w(l | 1123) = p 

w(l | 1122) =w{l | 1133) =p/2 

w(l | 2233) = (1 -p) 

w(l | 1222) = w(l | 1333) = w{\ \ 3222) = w(l | 3312) = (1 -p)/2 (8) 

The other rules are obtained by permutation of the neighboring sites and by 
cyclic permutation of the states. 

It is straightforward to check that the transition probabilities Wi(a), for 
both models, are invariant under the symmetry operations of the group C^ v . 

The nonequilibrium models, have the same symmetries as the Hamilto- 
nian of the three-state Potts model, given in (|4]) although in the present case 
the models are not defined by a Hamiltonian and do not satisfy detailed 
balance condition (|5|). That is, these models are microscopically irreversible. 
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3 Cumulants and scaling properties 

A convenient way to analyze the present models is through the use of the 
variables 

1 N 1 

X « = J^Y,( 6 ( a i> a ) - 3) ( 9 ) 

where a assumes the values 1, 2 and 3 and S(x,y) is the Kronecker delta. 
The order parameter has three components xi, x 2 , and £3 but just two of 
them are independent and the following property 

x 1 + x 2 + x 3 = 0, 

holds. 

It is useful also to introduce a set of homogeneous functions I n (xi, x 2 , £3), 
of a given order n, that are invariant under the symmetry operations R 
defined above. There is just one independent second order invariant given 
by 

I 2 = l(xl + x 2 2 + xl), (10) 
and just one independent third order invariant 

h = ~(xl + x 3 2 + xl), (11) 
The fourth order invariant function is 

h = l (xt + xi + xi) . (12) 
Again there is just one independent fourth order invariant. 

3.1 Cumulants 

In the present study the quantities of interest are the order parameter, 

m=\rn\ = (^y (13) 

and the reduced cumulant, 

^4 = ^ (14) 
6 



and, also, the Binder fourth order cumulant, which in terms of J4 and I2 
above defined is given by 

U=l- hj 2 ±. (15) 
We also analyzed the behavior of the order three reduced cumulant, 



= 777^ (16) 



3.2 Scaling properties 



Xo = — x y, 



The order parameter rrt has two independent components x and y and we 
will denote ^ 

= -y=(x~T + y~?) (17) 
where x and y are related to x±, X2 and 23 by the relations 

Zi = — Y X ~2 y ' ^ ' 

2 " 

x 3 = y. (19) 
The invariants fllPf ), (|ll|) and (|12"D are written, in terms of a; and ?/, as 

/ 2 = i(x 2 + y 2 ) (20) 

h = l(x*y- y j) (21) 

/ 4 = ^(x 4 + 2xV + Z/ 4 ) (22) 

Our main interest is to calculate the moments of the distribution asso- 
ciated to the order parameter. The moment M n of order n can be defined 

by 

M n =< \m\ n >= J \m\ n P(m,e, L)dm. 
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where P(m,e,L) is the probability distribution of m = \frt\, where fft = 
■j=(x,y), and e is equal to the deviation of the external parameter from its 
critical value and L is the system size. We assume that 

P(m, e, L) = L^ u (f)(m/e p , Le u ) . (23) 

Defining z = m/e 13 , 

M n , L = e P(n+l) L P/u J \ z \n^ ZiLe ^ dz _ 

Then, 

M n ,L = e^{e v Lf/ v F n {eL^) (24) 

where F n (Y) is an universal function. 

From relation (0), we get the following scaling forms 

M L =< m >= L- p/u m(tL l/u ) (25) 
A//" 4 ~ 

and 

/If 3 

±— = U 23 (eL^). (27) 



(Ml) 3 / 2 

where e is the deviation of the external parameter from its critical value 
and m(x), U 2 3(x) and U2i(x) are universal functions. 

For an infinite system the correlation length diverges as £ ~ e~ u and the 
scaling forms give the behavior m ~ e 13 for the order parameter. Moreover, 
the reduced cumulants U 2 i and t/23, defined in equations ( |X4|) and flTSD, are 



expected to attain according to (|26|) and ( p7|) a universal value at the critical 
point, that does not depend on the lattice size. The same behavior, of course, 
must hold for the reduced fourth-order cumulant U given in equation (|i~5|). 



4 Exact results 

When the temperature T — * 0, the equilibrium model defined by the Hamil- 
tonian (|J), will be in the ordered phase and with probability 1 in one of the 
three Potts states. For the nonequilibrium models, defined in (|7|) and (|8|), 
this limits correspond to p — >• 1. 
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In this limit it is expected the following behavior 



< M n >-> {(~x~r + (i)(-^r + 4(4 r} - (28) 

with a = 1, 2, and 3 and n = 1, 2, 3, Depending on the initial conditions 

the system will be with a* = 1, for all sites in the lattice, or in state where 
<7j = 2, for all sites in the lattice, or in the state o~j = 3, for all sites in the 
lattice. The factor (|) in equation takes into account this fact. For 
example, the second order invariant defined in expression (|T0|) will attain, 
the value 

J2 ^ l {( l )( l )2 + ( l )( ~l )2 + ( l )( ~l )2} 
that is J 2 -> 2/9 whenT -> (p -> 1). 

Following the same procedure, the limiting values of the third and fourth 
order invariants and reduced cumulants can be easily evaluated. In particular 
the values of the fourth order cumulant, defined in fll4|), in the limit T — ► 
(p — > 1), will be 

U 24 -+ 1. (29) 

which implies that the Binder fourth order cumulant, according to equation 
flT5|)) takes the limit, 

(3o) 

On the other hand when the temperature T — > oo (p — ► 1/3) the equi- 
librium system (the nonequilibrium symmetric model), is in the disordered 
state, 

Xi = X2 = £3 = 0. 

So the probability distribution associated to the order parameter is a distri- 
bution of independent variables, a Gaussian distribution, and we can write 

P(x,y) = ^exp[-(x 2 + y 2 )/2a], (31) 

with 

a=< (x 2 + y 2 )/2 > . 
So in this limit we have that the expressions (^Up and fl2"2p are related by 

h = 3(/ 2 ) 2 . (32) 
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Which implies that 

U 24 -> 3 (33) 

and 

U -»• 0. (34) 

The third order invariant /3 — ► when T —>■ oo (p — *■ 1/3). So C/23 — > 0. 

5 Monte Carlo simulations 

The system evolves in time according to the local rules and eventually reaches 
a steady state that can be of two types: a disordered steady state, where there 
is an equal average number of sites in each one of the three Potts states; either 
an ordered steady state characterized by the predominance of sites in one of 
the Potts states. 

The simulation of the equilibrium and the nonequilibrium models with 
C$ v symmetry was performed by considering square lattices with L 2 = N 
sites, and periodic boundary conditions. Each simulation started with a con- 
figuration generated at random and averages over several simulations where 
taken to get the final results. 



5.1 Equilibrium model 

We consider several values of the external parameter, the temperature T. We 
pick a site % at random and then we apply the Metropolis prescription to 
update site % according to the expression @ as follows. Let the state of site 
% be <7j. We change the site variable to Oi + a and calculate AH according 
to expression @ considering the nearest neighbor sites of the site % (which 
have not changed, since we are considering a one site dynamics). If AH < 0, 
then the new state will be a ia = (01, ...,ai +a , <x/v). Otherwise, if AH > 
0,we calculate p = exp(—(3AH) and generate a random number ( equally 
distributed in the interval [0,1]. If £ < p then the new state will be a 1 , 
otherwise the state does not change, that is, remains a. 

The system evolves in time until it reaches a stationary state. The time 
taken to the system to reach the stationary state depends on the temperature 
and on the lattice size. After discarding the first configurations, we used the 
following states in order to evaluate the state functions cumulants of the 
distribution probability associated to the order parameter. The stationary 
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states are equilibrium stationary states, i.e. they satisfy the detailed balance 
condition. As expected we found two types of stationary states: a ordered 
one where the order parameter m is different from zero and a disordered one 
where m — 0. 

The critical temperature for two dimensional three state Potts model 
is given by k B T c = l/(ln( v / 3 + 1) ~ 0.99497 @. In our simulation we 



take ks = 1 and we analyze the behavior of the cumulants (|14]) and fli~6|) 
as a function of the temperature and for different lattice sizes. As we can 
see in Figure 1 and Figure 2 when T — > T c the cumulant U24 — > U24 and 
the cumulant U23 — > C/| 3 , respectively, and it follows that t/| 4 and C/| 3 are 
universal. That is in the critical point U24 and U23 attain universal values that 
do not depend on the lattice size. We used the finite size scaling relations 



(|27| ) and (f26|). The expression (|27|) is related to the reduced Binder pO 
fourth order cumulant that so assumes an universal value U*. The value for 
these functions at the critical point where found by us to take 

C/ 2 * 4 = 1.16 ±0.01 (35) 

and {7* =0.61. The value of U23 at the critical point (see Figure 2) is U% 3 = 
0.245 ±0.01. 



5.2 Nonequilibrium models 

We consider several value of the parameter p. At each time step just one site 
is chosen at random and it is updated according to the prescriptions given in 
Section 2.2 are applied (rules (0) for the symmetric model and rules (§) for 
the majority model). After a transient, which depends on the model, on the 
size of the system and on the value of p, the system attains a steady state. 
Our simulations show that both models exhibit continuous phase transitions 
with the ordered steady state (m 7^ 0) occurring at high values of p. As p is 
decreased the transition takes place at a critical value p c , which is different 
for each model, and the system becomes disordered (m = 0) for p less than 

Pc- 

Using the finite size scaling relations (^7|) and (^) we obtain the criti- 
cal value p c for each model. For the symmetric lattice gas model, as shown 
in Figure 3 and Figure 4, the curves of f/ 24 versus p and the curves of U23 
versus p, for different values of L, intercept at the critical point p c estimated 
to be p c = 0.892 ± 0.003. It is worth to call the attention that this model 
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is similar to the one considered by one of us in a previous work [I7|. Both 
models evolves in time according to the same local Markovian rules. How- 
ever the one considered in reference [17] is a probabilistic cellular automaton 
(synchronous update) whereas the present model evolves in time according 
to a sequential dynamics (asynchronous update). As irreversible models are 
defined by the dynamics itself we do not have to expect the same value of 
the critical parameter. In fact they are different. We also observe that the 
results obtained previously |17j for the fourth order cumulant for the proba- 
bilistic cellular automaton are not sufficient precise due to presence of large 
fluctuations. Contrastingly, in the study of the stochastic lattice gas models 
considered here, we verified that the behavior of the cumulants, both of third 
and fourth order, are smooth. 

Figure 5 and Figure 6 show the curves of U24 versus p and the curves of 
U23 versus p for different values of L for the majority stochastic lattice gas 
model. The interception of these curves yields the critical point estimated 
to be p c = 0.883 ± 0.001. The values attained by U23 and U24 are universal 
and are f/ 2 * 3 = 0.244 ± 0.01 and U% 4 = 1.16 ± 0.01. And the Binder reduced 
fourth order cumulant is f/*=0.61. These universal values are in agreement 
with the results for the (equilibrium) three state Potts model. 



6 Summary 

We have considered systems that undergo a phase transition from a state 
with high symmetry to a state with lower symmetry. The phase with high 
symmetry is invariant under the symmetry operations of the symmetry group 
C% v . We analyzed first, the equilibrium three states Potts model, and then 
two irreversible models. All models have a continuous time evolution, gov- 
erned by a master equation. In the first case the model is defined by an 
Hamiltonian whereas in the second case they are defined only by transitions 
rates that do not obey the detailed balance condition. The phase transition 
that takes place, as an external parameter is varied, is a continuous phase 
transition from a disordered steady state to an ordered steady state. For the 
case of the equilibrium model, the phase transition occurs when the tem- 
perature T is varied and for the nonequilibrium cases when the parameter p 
is varied. We introduced a set of homogeneous functions that are invariant 
under the symmetry operations of the group 03^. From these functions we 
define the order parameter and the cumulants. The critical points were es- 
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timated by numerical simulations on regular square lattices of different sizes 
and using finite size scaling theory. Analyzing the cumulants we conclude 
that irreversibility plays an irrelevant role in the critical behavior and is not 
a property that might change the universal behavior. In fact, the fourth- 
order and the third-order cumulants attains universal values at the critical 
point and we found that this values are the same for the equilibrium system 
and for nonequilibrium systems, whenever periodic boundary conditions are 
considered. It is worthwhile to observe that the value of the cumulants may 
depend on the boundary conditions as was established for the equilibrium 
two dimensional models plf . 
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Figure Captions 

Figure l.The reduced cumulant U24 as a function of T for L = 14, 20, 28, 
and 40 (square lattices) for the three state (equilibrium) Potts model. 

Figure 2. The reduced cumulant U23 as a function of T for L = 14, 20, 28 
and 40 (square lattices) for the three state (equilibrium) Potts model. 

Figure 3. The reduced cumulant U 2 a as a function of p for L = 14, 20, 
28 and 40 (square lattices) for the three state symmetric model (nonequilib- 
rium). 
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Figure 4. The reduced cumulant U23 as a function of p for L = 14, 20, 28 
and 40 (square lattices) for the symmetric model (nonequilibrium). 

Figure 5. The reduced cumulant U 2 4 as a function of p for L = 14, 20, 28 
and 40 (square lattices) for the three state (nonequilibrium) majority model. 

Figure 6. The reduced cumulant U23 as a function of p for L = 14, 20, 28 
and 40 (square lattices) for the three state (nonequilibrium) majority model. 
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